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We consider the behavior of an Ising ferromagnet obeying the 
Glauber dynamics under the influence of a fast switching, random ex- 
ternal field. In Part I, we introduced a general formalism for describing 
such systems and presented the mean-field theory. In this article we 
derive results for the one dimensional case, which can be only partially 
solved. Monte Carlo simulations performed on a square lattice indicate 
that the main features of the mean field theory survive the presence of 
strong fluctuations. 

PACS numbers: 05.50+g 05.70.Jk 64.60Cn 68.35. Rh 75.10.H 82.20.M 

This paper considers the randomly driven Ising model (RDIM), introduced and 
discussed from a general point of view in Part la, in one and two dimensions. The 
main interest in studying such highly nonlinear nonequilibrium statistical physical 
systems lies in their possible applications for storage and other information pro- 
cessing tasks. Many biological systems, especially networks of real neurons, share 
common features with the RDIM in that they form strongly coupled systems driven 
by external stimuli acting in unison over a macroscopic number of elements in char- 
acteristic times shorter than the typical system relaxation time. Therefore, the 
Gaussian or Poissonian external "noise" is transformed radically within the coupled 
system, leading to a "correlated noise" with a lot of "strange" properties. The mean 
field theory was developed in Part la for a random binary switching external field. 
Applications to cortical neurons will be discussed elsewhere. 

The order parameter of the RDIM is a nonequilibrium stationary magnetization 
distribution, which undergoes a symmetry breaking bifurcation at some critical field 
strength. The analytic structure of this distribution also changes in character as 
a function of the temperature and field parameters. Hence, transitions between a 
singular function with fractal support to a singular function with euclidean support 
and further, to an absolutely continuous distribution can be observed. They can be 
seen in finite size effects and in the variance (the fluctuations) of the free energy. 
Similar transitions can be seen also in one dimension. However, the critical lower 
dimensionality of the RDIM remains one. Although the static critical exponent is 
a function of the field/temperature ratio, the dynamic exponent A = 2 remains 
unchanged. Some simple arguments will be given, supporting these results. 

The situation is much more difficult in two-dimensions, where only a Monte Carlo 
simulation approach was possible. The main difficulty is that computing the sta- 
tionary properties of the system requires a large population of different trajectories. 
We found that it is more convenient to use systems of moderate size than to rely 
on only few dynamic trajectories. Even so, the simulations require a vast amount 
of computer resources. We solved this problem by using the Siemens-Nixdorf neu- 
rocomputer SYNAPSE-1N110 for quite a different task than what it was conceived 
for. This machine is basically a matrix-computer, allowing us to run many dif- 
ferent systems in parallel. Even so, we were able only to determine roughly the 
phase diagram itself: we have no means at the moment for calculating the critical 
exponents. 

The paper is organized as following: in Section [| we consider the one dimensional 
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case. Although it cannot be fully solved, many interesting exact results can be 
derived. Section || deals with the Monte Carlo simulations we performed on finite 
square lattices. Many features of the mean-field dynamics are shown to survive the 
strong fluctuations characteristic to two dimensional systems. However, in contrast 
to the mean-field approach, the two dimensional model displays also an interesting 
spatial structure related to droplet dynamics. Comparisons between mean- field 
theory and two-dimensional results are systematically presented, including some 



preliminary results for hysteresis. Finally, we discuss our results in Section [II 



I. RDIM IN ONE DIMENSION 



A. The Master Equation 

In his pioneering work, Glauber! defined a stochastic dynamics where only single 
spin flips are allowed and hence neither the magnetization (the order parameter) 
nor the energy is preserved (model A universality class, seeu). Although intro- 
duced mainly for mathematical convenience, this dynamics is believed to describe 
appropriately many Ising-like systems. 

The energy of an Ising chain with periodic boundary conditions is given by 

N N 

i=i i=i 

Following n 

the notation of Part lH we denote by the vector jl = (si, S2, ■ ■ ■ , , — Sj, ■ ■ ■ , sn) a 
given configuration of spins and by /i, = (si, S2, • • ■ , — Sj, . . . , sjv) the same config- 
uration but with the spin Sj — > — Sj flipped. The external field is sampled from 
p(B) = \8(B — Bo) + \S(B + B ) at time intervals of length tb, 

00 

B(t) = p B Bp(B) 9(* " nT B)Q ((n + 1)tb - t) (2) 

n=0 

As explained in Part I, the Master Equation has the form 
P({ Si };t) = -C B(t) P({ Si };t) 

N N 

= ^2 w (P\^i) P (Pii t ) - BO 2 : f ) ( 3 ) 

i i 

where w(si) = w(jl\jli) is the transition rate from a configuration {si} into the state 
where only the i-th spin is flipped, {— Si}. 

Strictly speaking, the system we are going to consider here will never reach the 
equilibrium Boltzmann distribution e~@ E , where (3 — 1/ksT. Nevertheless, we 
require the detailed balance condition to be fulfilled for a constant value of the 
external field. 

The transition probability w(si) is determined by the constraint of detailed bal- 
ance only up to a positive arbitrary function of the neighboring spins, /(sj_i, Sj+x). 
Since in one dimension the phase transition is at T c = 0, the choice of the transition 
probability influences the analytic form of the critical singularitiesu. In what follows 
we will use the form 

w( Si ) = — [1 - Sl tanh(^ s i + H )] ( 4 ) 
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where K = (3 J, H = P^bB, a sets the time constant, and denotes next 

neighbor pairs. Glauber introduced this form for H = but used a slightly different 
one for H ^= 0. In one dimension one has then 

= tH 1 - s l ta,nh(K{s.^ 1 + s i+1 ) + H)] (5) 

For H = the Liouville operator Cb can be mapped onto a free-fermion spin- 
chain Hamiltoniam. This explains why the equations for the averaged spin prod- 
ucts, (IljeQ s j)' decouple in subspaces which can be diagonalized by appropriate 
Fourier transformations. As shown below, this property is inherited by the first 
moments, [7r Q (f)] of the stationary distribution V s . 



B. Magnetization and correlation functions 



Consider, for example, the time evolution of the local magnetization 

mi(t) = {si) t = J2 s i p ({si};t) (6) 

which can be obtained from (|^) as 
2 1 

= (siw(si)) = - [-mi(t) + (tanh(if ( Sl _! + s l+1 ) + H)) t ] (7) 



III: 



In order to make thfi-,relationship to chaotic maps more evident, we use now the 
'coarse grained' forrru of Formally, this procedure corresponds to a forward 
Euler discretization of ([?]), setting the time step equal to rg, and measuring time 
in units of a — tb- 

rrnit + 1) = (tanh(if + s i+1 ) + H)) t 

_mi-i(t)+m i+ i(t) 
= a + 7 ^ 1" o(sj-is i+ i) t (8) 

where 7 = tanh2if, h = t&nhH , and 

h 1 — 

°=2<1^ +1 ) (9) 
1-h 2 

^ = ^T-W (10) 

Similarly, for the correlation function 

= (siSj) (12) 

one obtains 

2 

Cj j = (s l s j (w(s i ) + w(sj))) t (13) 

a 

leading to the map 

Cij(t + 1) = i(s i tanh(i4T(s i _i + s l+1 ) + H) + s;tanh(iC(sj_i + Sj+x) + H)) t 

2 (14) 
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where the time unit is now set to a = 2tb- These recursions can be written again 
in terms of the variables a, b and 7, Eqs. (^|Tl|). p 

For H = one has 7 = 7 and, as shown by GlauberB, the slowest relaxation time 
equals the inverse of the smallest eigenvalue of the magnetization subspace, Eq. (Q). 
This relaxation time diverges with the square of the static correlation length (the 
dynamic critical exponent is z = 2). 

For the binary field distribution of (Q) the map of the local magnetization has 
two branches: 



mi(t + l) 



a+ ~ m i - 1 (t)+m i+1 (t) +6c ._ M+1 ft) withprob. i 

- fl + 7 m - l(t) + mi+l(t) -6 Cl _ M+1 (f) withprob. i (15) 



where it is implicitly assumed that the time needed to switch the field is negligible, 

Tswitch «Tb- 

From Eq. ([l^) it is evident that the map for the local magnetization couples to a 
two-spin correlation, which in turn couples to higher order correlations, eta. Hence, 
the full dynamic map lives in a 2 N dimensional space, as stated in Part iEl. 

Nevertheless, some partial results can be obtained for the stationary distribution. 
Define a 'thermal' and a 'dynamical' average, (...) resp. [...]. In the stationary 
state one has [(A(t + 1))] = [(A(£))] for any spin-function A. The average of the 
local magnetization obeys 

[7] 

\ m i\ = "y(L m i-l] + [ m i+l]) (I- 6 ) 

where (recall that 7 = tanh2i4T, ho = tanhiJo) 

H = ^W (17) 

For the translation invariant magnetization m = ^ i rrii and the two-spin cor- 
relation function c m = n jf 53j c M+m on e obtains equations formally similar to the 
ones solved by Glaubera. From (|l6|), except for [7] = 1, the magnetization vanishes. 
The stationary two-spin correlation function obeys 

[Cm] = l(W + [<Wl]) ( 18 ) 
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which leads to 



[c ro ]=r ? H with v = L_^L_J2l! (19) 

m 

In general, when expressing the Frobenius-Perron operator in the basis formed 
by all moments of spin-products [7f ■?(*)], the subspace of the first moments (q = 1) 
is closed and can be diagonalized by a Fourier transform. The remaining part, 
however, is intractable. For example, the second moment of the magnetization 
reads 

[to 2 ] = [a 2 ] + 2 [aft] [c 2 ] + [ 7 2 ] [to 2 ] + [6 2 ] [c 2 ] (20) 

and couples both to the first and to the second moment of the translation invariant 
correlation functions, [02] and [c 2 ,], respectively. 
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C. The T c = phase transition 



As expected, in the stationary state the odd moments of the magnetization vanish 
at T > 0. By expanding 7 and 77 at low temperature one obtains after straightfor- 
ward calculations that close to T c = the correlation length is in leading order 

if 2K > H 

(21) 

if 2K < H 

where Hq = PpbBq and Bq is defined as in Eq. (||). Hence, in one dimension the 
RDIM has a critical line with continuously changing singularities at T c = for 
k = Hq/2K G [0, 1]. The slowest relaxation time corresponds to the magnetization 
(order parameter) decay and can be computed as 

( t 2 if k < 1 
r- y l = 2(1 - [7]) ~ { (22) 
[2 if k > 1 

Consider first the case of a strong field, k > 1. The field will align all spins in one 
iteration step, as evident from Eq. (p2[). Hence, the spins are almost always parallel 
to the driving field. Since [7] vanishes, there is no spontaneously broken symmetry 
and [to] = due to the symmetry of the field distribution, p(B) = p{—B). 

For fields smaller than the critical field K c = 1 one obtains a true symmetry 
breaking ferromagnetic phase. Interesting enough, while the divergence of the cor- 
relation length decreases continuously according to Eq. (|2l]), the critical dynamic 
exponent remains z = 2 up to and including k < 1. A physical argument explaining 
this result is presented below. 



In 77 
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D. Kink dynamics at T = 

Consider now the transition probability w(si) at T = 0, Eq. (^), which is a 
function of the three spins Sj, Sj-i, and Sj+i. Let us call an interface between two 
oppositely oriented spin domains a kink. If k — Hq/2K = and T = each kink 
performs a random walk, moving with equal probability to the left or to the right. 
When two kinks become nearest neighbors, they annihilate because in the next time 
step the single spin left between them will flip with probability one. 

Due to this annihilation process the number of kinks decreases steadily and in 
the end only very few are left. Two kinks situated at the typical distance £ (the 
correlation length) will meet via diffusive motion in the characteristic time r ~ £ 2 , 
which explains why the critical dynamic index is z = 2. The situation is similar if 
one is close to (but not at) T c = 0. 

How does this picture change if we switch on the random external field [k > 0)? 
The domains parallel to the external field start growing - both ends of such a cluster 
will move outwards. Once the field changes sign, these domains shrink again and 
the clusters of oppositely aligned spins grow. During a longer period of time these 
effects compensate each other and the surviving kinks perform effectively a random 
walk. However, the annihilation rate of kink-pairs is highly increased. Assume, for 
example, that after 2T iteration steps the external field had the value +Hq (T + n)- 
times. The probability for this to happen is given by the Bernoulli distribution, 
Bn T — CrK^)™' D ur hig this time, a down-oriented domain whose original length 
was Lq shrinks on average by 2n(2p— 1), where both kinks associated with the ends 
of the domain have moved inward during each of the n time steps with probability 



5 



p > 1/2 and outward with probability q = 1 — p. Hence, if 2n(2p — 1) ~ Lo, the 
domain will be eliminated, p is a function of re = |j8>, e.g. p(re = 0.2) = 0.982. 
Due to the depletion of small clusters, the number of kinks decreases much faster 
than in the absence of the external field. This effect is illustrated by a numerical 
simulation in Fig. [jj showing the T = dynamics of walls for k — and k = 0.2, 
respectively. 

Once only a few large clusters remain, however, their width becomes macroscopic. 
On this scale the kinks perform again a random walk and asymptotically one regains 
z = 2. 

At finite temperatures, however, the presence of the field term facilitates the 
nucleation of new clusters, so that the correlation length ( pl| ) (the mean cluster 
size) is less divergent when Hq > 0. 

E. The magnetization distribution 

Consider again the map (|l5|). In the translational invariant sector one has 

!a + jm(t) + bc%(t) with prob. \ 
(23) 
—a + ym(t) — bc2(t) with prob. ^ 

As already discussed, the magnetization couples to the correlation function C2 (t) , 
etc. A simple approximation to decouple the magnetization sector is using for C2(tX 
the stationary value C2 = rj 2 . The resulting map corresponds to a Bernoulli-shiftfl 
and is shown graphically in Fig. ||. If the gap between the two branches is positive 
A > 



A= 2( fl + 6c a )(l-27) >0 

1-7 

the corresponding stationary magnetization distribution is a Cantor set. 

Since a, 6, (1 — 7) are positive and ci > for ferromagnetic interactions, the 
demarcation line between a fractal and a nonfractal magnetization distribution is 
given by 

1=\, (25) 

independently of the actual value c^if) might have. The time dependence of C2(t) 
induces nonlinearities in the map. Therefore, although the distribution remains 
fractal for 7 < 1/2, in general it is not a homogeneous Cantor set. 



II. MONTE CARLO SIMULATIONS IN TWO DIMENSIONS 



We simulated the RDIM on a two dimensional square lattice on the neurocom- 
puter SYNAPSE-1/N110. In the following sections we first describe a Monte Carlo 
Algorithm (MCA) designed to make use of the computational power of SYNAPSE- 
1 and then present our numerical results. They include a phase diagram in the 
H-K-plane, magnetization distributions in the para- and ferromagnetic regime, and 
a series of snapshots documenting the behavior of the system. 
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A. The algorithm 



SYNAPSE-1 is a workstation-driven coprocessor consisting of a systolic array 
of eight MA16 Neural Signal Processors. It was kindly put at our disposal by 
the ZFE of Siemens AG. Its hardware was designed to tackle typical problems 
encountered when simulating neural networks, namely calculations involving very 
large matrices. In order to efficiently make use of the CH — \~ library interface supplied 
with SYNAPSE-1 for the RDIM, we devised a MCA that simulates multiples of eight 
lattices in parallel. 

Consider a square lattice of spins g {—1,1} of linear dimension L, where k 
numbers the system and i,j — 1 . . . L denotes the lattice position. Each of eight 
lattices is represented as a column vector s v by renumbering indices v = iL + j. 
The eight systems can thus be treated as one 8 x L 2 -matrix. By setting v i— > v + L 2 
for v < and v i— » v — L 2 for v > L 2 we enforce helical boundary conditions. The 
neighbors of spin v are [i — v ± I, /i = v ± L. 

In order to avoid metastable states induced by a simultaneous update of neighbor- 
ing spins we split the lattices into black and white sites in a checkerboard fashion, 
leading to two matrices encoding the eight systems. Note that when a lattice is 
divided up in this fashion, if L is odd, the sites in the first and last row have neigh- 
bors of their own color. If L is even, the same is true for the first and last column 
of the lattice. For technical reasons we chose L to be odd. 

To further simplify the updating scheme, a copy of the first and last L components 
of each lattice are included at the end respectively beginning of each column vector. 
A Monte Carlo step now consists of a parallel update of all black sites followed by 
an update of all white sites (or vice versa). 

From the well known Glauber dynamic rule Eq. (Q), a spin is flipped under the 
following condition. Given a random number z S [0, 1] drawn from a uniform 
distribution, spin s„ is updated according to 

( -s u if z< |(l-tanh(s I/ (JCS (tl s /J + if))) 
s v ^> I ' (26) 

[ s v if z > i(l - tanh(s w (ifS (U s M + H))) 

Usually, spins are either treated sequentially or chosen randomly. SYNAPSE-1 
permits the parallel generation of (pseudo-) random numbers in a single Elementary 
Operation (ELOP) which can be piped through a function lookup table at no extra 
computational cost. For this reason we transform the flip condition Eq. ( p6|) into 

!-s v if EpS^ < 2^ log(|^7) ~ s » 7c 
(27) 
s v if S^s^ > lo S(f^7) - 

where z' is drawn from a uniform distribution in [— ^, i]. The RHS of the flip 
condition is evaluated in two ELOPs: One to generate a matrix of random numbers 

piped through the function z' i— > ^ log( f— p ) and one weighted matrix subtraction. 

The LHS also takes two ELOPs to calculate from the "black" and "white" matrices. 
Two further ELOPs are required to construct and evaluate a flip indicator matrix. 
Some more operations are necessary to fix boundary conditions and to evaluate 
the mean lattice magnetizations. This procedure is applied sequentially first to 
the matrix holding the "black" spins and then to the "white" one to accomplish a 
complete Monte Carlo step. 

Initially, the spins in a lattice are set to +1 with probability p and —1 with 
probability 1 — p, where different values of p can be used for each lattice in one 



7 



simulation run. The results for systems of linear dimension L = 415 are initialized 
withp =0, 0.2, 0.4, 0.5, 0.6, 0.8, and 1, in addition, there is one lattice in which the 
top half of all spins is set to +1 and the bottom half to —1. The external driving 
field is the same for each system. Simulations at smaller L consist of 64 systems 
initialized with p — \, but each with its own driving field trajectory. 

Temperature K is measured in units of the critical temperature, K — > where 

K c = -^jr = \ h^v^ + 1) « 0.44069 of the standard two dimensional Ising model, 
i.e. K — 1 corresponds to the critical temperature for Ho = 0. 

B. Dynamics and phase diagram 

In order to understand the dynamics of the two-dimensional RDIM in more detail, 
it is useful to consider first what happens to a cluster of parallel spins at T = 0, 
in analogy to kink dynamics for the one dimensional case. Transition rates at 
T = are either 0,1, or |. Recalling that K is measured in units of K Cl define 
k' := = kK c . Consider now a square cluster of 2N x 2N parallel spins under 
the influence of an anti-parallel external field. First, if < k' < 1, the spins at the 
corners of the cluster flip with probability p — 1 , all other remain anti-parallel to the 
field (as shown in Fig. [| a) . Thus the cluster disappears if the external field remains 
constant for 2N — 1 consecutive steps. Secondly, if 1 < k' < 2, such a cluster will 
be destroyed in N steps, due to the fact that all but inner spins will flip with p = 1. 
Thirdly, at k' > 2 we arrive at a driven paramagnetic phase. Regardless of their 
position, all spins will flip into the direction of the driving field with p — 1. For 
the case of k' = 0,k' = l,or k' = 2, corner, edge, and inner spins flip with p = i 
This implies that, e.g. for a strong driving field with k' = 2K C , nucleation flips (see 
Fig. |3| c) may take place inside the cluster, creating magnetic swiss cheese. 

The transition probabilities of the processes shown in Fig. || increase with increas- 
ing field strength. At small fields the a)-type flip is prevalent, resulting in a radial 
growth (shrinkage) of clusters. Although energetically more expensive, the b)-type 
flip has a large entropy contribution and results in long-wavelength growth of flat 
domain walls. The nucleation process shown in flip c) has the smallest probability. 

What happens when switching instantly the field from the equilibrium state 
at — Hq into the unfavorable direction Hq? The system relaxes from the now 
metastable state to the new equilibrium value. Obviously, the lifetime of the 
metastable state depends on the strength of the applied field. This scenario has 
been discussed in detail in the ferromagnetic phase using droplet theory and Monte 
Carlo simulations (see, for examplejj). Here, four distinct field intervals, shown 
schematically in Fig. ^ were, identified in which the lifetimes markedly differ due 
to different decay mechanisms. A numerical result for K = 1.25 is shown in Fig. 
H where we approximated the metastable lifetimes by measuring the average first 
passage times (FPT) from m = —1 to m = 0.7 in Monte Carlo steps. Figs, g and ^| 
show examples of the time development of the magnetization from the metastable 
state to th new equilibrium for the mean field and 2D model. 

We calculated numerically a phase diagram in the K — j^— plane for the 2D 
RDIM, Fig ||. Similar to the mean field and one dimensional model, there is a 
paramagnetic, a ferromagnetic (and a driven paramagnetic) phase. Note that for 
K i — > oo, the phase boundary should remain below k! = 2. We are currently in no 
position to assess this. Also, the behavior around K = 1 doesn't seem to correspond 
to the first order dynamic freezing transition seen in the mean field theory, rather 
a second order transition is likely. 
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C. The paramagnetic phase 

If the external field is above its critical value, Ho > H c , the stationary phase of 
the RDIM is paramagnetic. In this phase the system relaxes relatively fast to the 
equilibrium state, except close to H c , where critical slowing down sets in due to the 
type-I intermittency effects discussed in the previous section. 

This behavior of the average magnetization is shown in Fig. and is probably 
enhanced by local correlations not taken into account in the mean field approxi- 
mation. If the random field is switched on, close but above H c the critical slowing 
down is dramatically enhanced. Further away from the phase transition point the 
dynamics is - similarly to the one-dimensional case - determined by the nucleation 
and radial growth (shrinking) of droplet-like clusters. 

As expected from the mean-field results, the RDIM can display a fractal magne- 
tization distribution at higher fields. This is shown in Fig. |^. Thermal fluctuations 
and finite size effects wash out the fine structure of the multifractal magnetiza- 
tion distribution predicted by the mean-field calculations. However, the presence 
of sharp peaks in the distribution (and their scaling behavior) demonstrates that 
some of the main features of the magnetization distribution survive the thermal 
fluctuations. 

These peaks are related to long-lived droplets whose radius is large enough to 
allow them to stay alive even when a long series of unfavorable external field draws 
makes them shrink. It is, however, the competition between the two thermodynam- 
ically stable states which leads to a chaotic dynamics and strange attractors. 

D. The ferromagnetic phase 

The situation is even more complex below H c . The schematic dependence of the 
average time spent in a thermodynamically unstable state vs. the inverse of the 
field strength is shown aftenj in Fig. [|. 

The external-field sampling time tb should be chosen in either the strong or multi- 
droplet regime (certainly not in the coexistence regime). By varying tb is seems 
possible to explore these different dynamic mechanisms in more detail. Fig. [To| 
shows a simulation in the ferromagnetic phase where one can observe both a multi- 
droplet (series 1 and 3) and a domain-wall type (series 2) dynamics. Again, the 
spontaneous magnetization distribution shows well separated peaks, which can be 
seen in Fig. [H]. 

For the sake of completeness, we show also a Monte Carlo simulation of the 
hysteresis measurement described in Fig. 8 of Part Sa. Only the evolution of one 
system initialized with spins up/down with equal probability is displayed. Again, 
one can see that the thermal fluctuations are smoothing only the fine scale structure 
of the mean- field predictions - but the main features remain intact. 

The results presented here leave open many questions regarding the 2D RDIM - 
an in-depth study by Monte Carlo simulation lies currently beyond our means. The 
theory and evolution of droplets and domains in the RDIM as well as hysteretic 
effects remain interesting research topics. 

III. SUMMARY AND DISCUSSION 

In Part I and in this article we have discussed in detail the properties of a simple 
Ising model subject to a fast switching external field. In many ways, the situation 
is just the opposite as in quenched random systems. While there the defects and 
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hence the (local) fields are frozen relative to the spin degrees of freedom, which 
are (in principle) free to relax, in the RDIM the external field is the fast variable 
compared to the interacting spin system. From a "technological" theoretical point 
of view the situation is, however, much easier. We expect that similar analytic 
results can be obtained for other random distributions as well. Strongly driven 
systems show rather peculiar properties, which one could use for increasing the 
storage properties of ferromagnetic materials. For example, it seems possible to 
use arithmetic coding in 'preparing' a semi-macroscopic ferromagnetic region to fall 
into a given distribution peak, as the ones shown in Fig. |p. An appropriately 
sensitive reading head could then discriminate between the different magnetization 
values. Such 'devices' could be tested first with the help of Monte Carlo simulations. 
However, the main application domain for randomly driven systems might well be 
in biology. Further work in that direction is under progress. 
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FIG. 2. The truncated ID magnetization map for A > 
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FIG. 3. The basic spin-flip configurations, Ho points upwards, up (down) spins are 
black (white), a spin turning up is drawn in grey, a) a 'droplet '-flip, b) a 'domain'- flip, c) 
nucleation-flip 
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FIG. 4. The average time spent in the metastable thermodynamic state as a function of 
inverse field-strength. The different domains are denoted according to their main relaxation 
mechanism. 
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FIG. 5. The average first-passage time from m — — 1 to m — 0.7 at K — 1.25 as a 
function of inverse field-strength. < t(m — 0.7) > is calculated from an ensemble of 64 
systems with linear dimension L = 143. 
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FIG. 6. Mean field iteration m vs. time t in an unfavorable field for = 0.057 and 
K = 1.2, which is slightly above the critical field. 
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FIG. 7. Monte Carlo simulation in an unfavorable field for = 0.08, K = 1.2, tb = 1, 
and linear dimension L — 143. The solid line shows the average magnetization in an 
ensemble of 64 systems, m, as a function of time t. The dashed lines mark a 1-a range 
around the mean. All systems start with all spins down, i.e. Vu : s v = — 1. 
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FIG. 8. Phase diagram of the 2D RDIM from Monte Carlo simulation of 64 systems of 
linear dimension L = 63. 
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FIG. 9. Top: Magnetization distribution for the square lattice RDIM averaged from eight 
different initial conditions for K = 0.4, -2. = 0.5, tb = 1, and linear dimension L — 415. 
Here, the simution covers more than 2 • 10 5 Monte Carlo sweeps. Note the similarity to 
Fig. 1 of Part I. Bottom: K = 1, ^ = 1, tb = 1, and L = 415. The simulation run covers 
close to 2 • 10 5 Monte Carlo sweeps. Compare this to Fig. 2 of Part I. 
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FIG. 10. Time series of the mean magnetization of the RDIM. Parameters are K = 2, 
= 1, tb = 1, and the linear dimension is 128. Note that here the boundary conditions 
are periodic. Snapshots of the system are shown in series 1, 2, and 3 on the next pages. 
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Series 2, continued 
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FIG. 11. Top: Magnetization distribution for the square lattice RDIM averaged from 
eight different initial conditions for K = 2, ^ = 1, tb = 1, and linear dimension L — 415. 
We simulated more than 2 • 10 5 sweeps. Only the region from m = 0.8 to m = 1.0 is shown, 
the distribution is symmetric in m. It is similar to the mean field distribution around the 
critical driving field H c (c.f. Fig. 3 and 4 of Part I). Bottom: K = 1.33, = 0.75, tb = 1, 
and L — 415. We tracked 1.5 • 10 5 Monte Carlo sweeps. Here, the distribution resembles 
the mean field distribution close to the critical driving field He- 
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FIG. 12. Hysteresis in the square lattice RDIM, showing the mean magnetization m(t) 
vs. the external driving field H(t), see Eq. (32) of Part I. Parameters are ^ = = 1, 
K = 2, = y§§q, and linear dimension L = 415. The simulation ran for more than 10 5 
sweeps. See Fig. 8 of Part I for comparison. 
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